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ABSTRACT 

The  experimental  carbon  Auger  lineshape  for  graphite  has  been 
obtained,  corrected  for  the  effects  of  the  secondary-electron 
background  and  extrinsic  losses  and  placed  on  an  absolute  energy 
scale  through  the  use  of  photoelectron  measurements.  The  resulting 
lineshape  is  compared  to  a  model  which  consists  of  the  self¬ 
convolution  of  the  graphite  one-electron  density  of  states 
including  atomic  values  for  the  symmetry-determined  Auger  matrix 
elements.  A  poor  comparison  results  from  this  simple  model  which 
is  considerably  improved  by  the  inclusion  of  dynamic  initial-state 
screening  effects.  Further  improvement  results  from  accounting  for 
final-state  hole-hole  interactions.  The  final  state  is 
characterized  by  effective  hole-hole  Interaction  energies  of  2.2 
eV,  corresponding  to  two  holes  in  the  o  band,  1.5  eV  for  one  hole 
in  the  o  and  one  in  the  u  band,  and  0.6  eV  for  both  holes  in  the  it 
band.  The  remaining  discrepancies  in  our  model  comparison  are 
suggested  to  be  due  to  a  plasmon  emission  intrinsically  coupled  to 
the  Auger  final  state. 
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I.  INTRODUCTION 


The  use  of  detailed  Auger  spectral  lineshape  analysis  to 

W 

obtain  local  electronic  structure  information  has  had  increased 
emphasis  over  the  past  few  years,  as  is  evident  from  the  number  of 
recent  review  articles  devoted  to  this  subject  {1-20}.  This 
interest  stems  from  the  local  nature  of  the  Auger  process  which  has 
as  its  initial  state  a  missing  core  electron.  For  core-valence- 
valence  Auger  transitions,  the  core  hole  state  captures  a  valence 
electron  and  transfers  its  excess  energy  to  the  ejection  of  a 
second  valence  electron,  the  measured  Auger  electron.  The  kinetic 
energy  (KE)  of  the  ejected  Auger  electron  can  be  approximated 
18,9},  by  the  expression 

KE  -  Ic  -  Ij  -  ‘k  -  Uerf.  (I) 

where  the  I's  are  the  one-electron  binding  energies  of  the  core  (c) 
and  valence  (J.k)  states  involved,  and  Ueff  takes  into  account  the 

interaction  between  the  two  final-state  holes.  Equation  1  of 
course,  refers  to  a  single  Auger  transition  while  the  Auger 
spectrum  is  composed  of  all  possible  Ij,  Ik  combinations.  This 

procedure  amounts  to  taking  the  self-convolution  of  the  set  of 
valence  states  Ij  or  1^,  in  other  words,  to  a  self-convolution  of 

the. density  of  states  (DOS).  The  local  nature  of  this  process 
stems  from  the  limited  spatial  extent  of  the  core  wave  function 
which  assures  that  the  Auger  process  probes  the  valence  electron 
density  over  the  same  spatial  extent.  The  implications  of  this 
local  sensitivity  with  respect  to  molecules  have  been  developed  in 
a  recent  review  {10}. 


The  C(KVV)  Uneshape  of  graphite  (the  notation  KVV  indicates 
that  the  core  hole  is  in  the  K  level  and  both  final-state  holes  are 
in  valence  levels)  has  been  the  subject  of  considerable  recent 
study  {21-261.  Although  much  of  this  attention  has  been  in  the 
context  of  studying  the  more  novel  graphite  intercalation  compounds 
{21,22,27-291  the  graphite  Auger  spectrum  is  itself  of  interest 
since  it  represents  the  infinite  limit  of  the  fused  ring  series: 
benzene,  naphthalene,  phenanthracene ,  etc.  In  this  role  the  C(KVV) 
lineshape  of  graphite  i3  unique  among  the  ring  aromatic  Auger 
llneshapes  because  the  two  final-state  holes  resulting  from  the 
Auger  process  have  a  chance  to  delocalize  over  a  much  larger  volume 
than  would  be  permitted  by  the  finite  size  of  the  molecules  {10). 
Thus,  it  is  possible  that  final-state  hole-hole  correlation  effects 
may  be  negligible  if  the  holes  actually  are  able  to  delocalize. 

In  addition,  graphite  is  a  model  system  for  studying  initial-state, 
core-hole  screening  effects  in  aromatic  systems.  Previous 
theoretical  calculations  have  indicated  that  core-hole  screening 
significantly  alters  the  shape  and  magnitude  of  the  measured  it  DOS 
{26,30),  but  the  effects  of  these  changes  in  the  graphite  Auger 
lineshape  have  not  been  examined. 

The  first  attempt  at  obtaining  an  accurate  C(KVV)  lineshape 
for, graphite  was  reported  by  Smith  and  Levinson  {231.  They 
utilized  a  data  reduction  procedure  which  has  become  almost  the 
standard  treatment  of  Auger  data  (31 1  in  order  to  obtain  detailed 
electronic  information.  The  data  was  taken  in  the  derivative  mode 
and  numerically  integrated.  A  background  was  removed  in  a  manner 
developed  by  Sickafus  {32-3*0,  and  the  resulting  Auger  lineshape 
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was  loss  deconvoluted  utilizing  a  263  eV  electron  elastic  peak  and 
attendant  loss  spectrum. 

An  attempt  at  quantitatively  interpreting  the  C(KVV)  lineshape 
for  graphite  reported  by  Smith  and  Levinson  {23)  has  recently  Keen 
reported  by  Murday,  et  al .  {21).  They  deduced  the  one-electron 
partial  DOS  (o  ,o  ,it  )  for  graphite  from  X-ray  emission  spectra 

3  p  p 

(XES),  X-ray  photoemission  spectra  (XPS),  and  an  assumed  electron 

2 

configuration  of  sp  ir.  The  Auger  lineshape  was  then  produced  from 

a  fold  of  these  one-electron  partial  DOS  assuming  noninteracting 

final-state  holes  and  no  screening  effects.  However,  an  error  (to 

be  discussed  later)  in  their  self-fold  makes  the  agreement  with 

Smith  and  Levinson  fortuitous. 

We  have  obtained  graphite  C(KVV)  spectra  which  show 

significant  differences  from  that  reported  by  Smith  and  Levinson 

{231  and  demonstrate  that  these  differences  are  due  to  an  improper 

loss  deconvolution  of  their  experimental  data.  Thi3  improper  data 

handling  resulted  in  incorrect  assumptions  in  the  subsequent 

theoretical  analysis  of  Murday,  et  al .  (21). 

Our  C(KVV)  lineshape  for  graphite  was  corrected  for  both  the 

effects  of  the  secondary-electron  background  and  the  extrinsic 

losses  suffered  by  the  Auger  electrons  in  leaving  the  solid. 

% 

Extrinsic  losses  are  those  external  to  the  Auger  process  such  as 
those  that  result  from  an  electron  moving  through  a  solid.  In 
contract,  intrinsic  losses  are  associated  with  the  Auger 
transition.  The  raw  Auger  data  were  taken  in  two  separate 
laboratories  and  on  three  distinct  types  of  electron  energy 
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analyzers  and  data  reduction  procedures  were  applied  independently 
In  each  case.  The  absolute  energy  scale  for  the  Auger  line  has 
been  established  by  utilizing  the  valence,' core  level  and  Auger 
features  available  concurrently  in  the  graphite  XPS  spectrum.  The 
resulting  llneshape  is  compared  to  a  model  which  consists  of  the 
self-convolution  of  the  graphite  one-electron  density  of  states 
Including  atomic  values  for  the  symmetry-determined  Auger  matrix 
elements.  A  poor  comparison  results  from  this  model  which  is 
considerably  improved  by  the  Inclusion  of  dynamic  screening  effects 
in  the  initial-state.  Further  Improvement  in  the  model  results 
from  accounting  for  the  effect  of  final-state  hole-hole 
Interactions  through  the  use  of  a  formalism  developed  by  Cinl 
(35,361  and  Sawatzky  (371.  From  the  resulting  parameterized-model 
fit  to  the  experimental  llneshape,  we  find  effective  hole-hole 
interaction  energies  of  2.2  eV  for  two  holes  in  the  0  band,  1.5  eV 
for  one  hole  in  the  0  and  one  in  the  it  band  and  0.6  eV  for  both 
holes  in  the  it  band.  Areas  of  discrepancy  remain  after  the 
inclusion  of  corrections  for  these  two  effects  of  many-electron 
processes  and  ' >ese  descrepancles  are  suggested  to  be  due  to  a 
jjiusinc..  .,«..Jon  intrinsically  coupled  to  the  Auger  final  state. 

II.  EXPERIMENTAL 

,  Evaluation  of  the  quality  of  the  data  presented  here  with 
respect  to  instrumental  and  possible  sample  differences  has 
Involved  a  "round  robin"  analysis.  The  C(KVV)  and  loss  spectra 
were  obtained  In  the  N( E)  mode  at  two  different  laboratories  on 
three  different  types  of  electron-energy  analyzers:  a  narrow- 
aperture,  retarding,  partial-hemispherical  sector  instrument;  a 
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medium-aperture,  retarding,  full-hemispherical,  sector  Instrument; 
and  a  retarding,  double-pass,  cylindrical-mirror  analyzer. 

Spectra  were  obtained  from  the  basal  plane  of  single  crystal 
graphite  (SCG),  highly  oriented  pyrolytic  graphite  (HOPG)  and  POCO 
graphite  (a  machinable  amorphous  graphite).  For  the  HOPG  and  SCG 
samples  a  layer  was  peeled  off  just  prior  to  insertion  i: to  the 
vacuum  system  and  all  samples  were  heated  to  -1200  K  by  electron 
bombardment  from  the  backside.  In  the  particular  case  of  the  SCG 
sample,  XPS  analysis  yielded  no  detectable  core-level  peaks  other 
than  C(1s)  even  after  several  days  in  the  vacuum  system.  The 
medium-aperture,  retarding,  full-hemispherical,  sector  analyzer 
system  had  an  associated  X-ray  source  and  a  Helium  lamp  In  addition 
to  the  electron  gun.  As  a  result  it  was  possible  to  obtain  in  the 
same  scan  the  X-ray  excited  Auger,  core  level,  and  valence  spectra 
all  calibrated  to  the  graphite  Fermi  level;  and  to  obtain,  without 
repositioning  of  the  sample,  the  ultraviolet  valency  spectrum  or 
the  electron-excited  Auger  spectrum. 

III.  RESULTS  AND  ANALYSIS 

Figure  1  shows  an  example  of  the  electron-excited  Auger 
results  from  the  POCO  graphite  sample  prepared  in  the  manner  just 
described  taken  In  the  N(E)  mode  with  the  narrow-aperture  analyzer. 
Raw  data  of  the  type  shown  in  Fig.  la  results  from  the  true  Auger 
spectrum  as  well  as  contributions  from  two  factors:  1)  the  spectra 
reside  on  the  secondary  electron  background  and  2)  the  outgoing 
Auger  electrons  lose  energy  in  traveling  through  the  near  surface 
region  of  the  sample.  The  second  of  these  effects  is  evident  in 
the  structure  and  the  featureless  tail  on  the  low  kinetic  energy 
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aide  of  the  spectrum.  To  obtain  a  more  representative  measure  of 
the  undistorted  Auger  lineshape,  these  effects  must  be  removed  from 
the  raw  data.  Houston  (31)  has  treated  these  lineshape  distortions 
in  detail  in  a  recent  publication  in  which  data  similar  to  that  in 
Fig.  la  was  used  as  an  example.  Various  analytic  functions  are 
described  for  direct  background  subtraction  and  the  effect  of  the 
losses  is  removed  by  a  deconvolution  procedure  utilizing  model  loss 
functions . 

In  applying  these  procedures,  it  is  assumed  that  in  the 
absence  of  the  background  there  exists  a  "true"  Auger  lineshape 
covering  a  limited  energy  range,  i.e.,  an  Auger  line,  which  when 
convoluted  by  the  proper  loss  function  for  the  material  under  study 
yields  the  background-corrected  raw  data.  The  loss  function  is  the 
lineshape  that  would  be  measured  if  one  had  a  monoenergetic , 
internal  source  of  Auger  electrons.  Such  a  source  is,  of  course, 
not  available  and  one  must  use  approximate  loss  functions  to 
accomplish  the  lineshape  correction.  An  approximate  loss  function 
that  is  easily  obtained  and  often  used  is  that  of  the  elastic  and 
near  elastic  electron  backscatter  spectrum  resulting  from  a 
monoenergetic  electron  beam  incident  on  the  sample  surface  at  an 
energy  near  that  of  the  Auger  feature.  Although  this  model  loss 
function  may  in  some  cases  not  be  adequate,  we  justify  its  use  here 
because  of  the  similarity  in  shape  that  we  find  between  the 
backscatter  spectrum  and  the  XPS  core-level  spectrum  for  graphite-- 
the  photoelectron  being  a  much  better  approximation  to  an  Auger 
electron  than  a  backscatter  electron.  The  principal  remaining 
difficulty  in  using  the  backscatter  approximation  for  the  model 


1 033  function  is  that  we  have  no  valid  model  for  the  "intrinsic 


losses  associated  with  the  Auger  process  and  these,  if  appreciable, 
remain  in  the  deconvoluted  Auger  lineshape.  This  point  will  turn 
out  to  be  important  in  the  ensuing  discussion. 

The  backscatter  approximation  to  the  loss  function  i3  show.-,  in 
Fig.  1e  taken  at  an  incident  electron  energy  of  290  eV.  One  can 
see  from  this  data  the  two  graphite  characteristic  losses  (usually 
described  as  "plasmon"  losses)  at  loss  energies  of  approximately  7 
and  27  eV;  these  appear  very  similar  to  the  previously  published 
(38)  loss  features  associated  with  the  XPS  C(1s)  line.  It  should 
be  noted  that  the  instrument  response  of  the  electron  energy 
analyzer  used  to  obtain  the  data  of  Fig.  1  is  given  by  the  shape  of 
the  elastic  portion  of  this  spectrum. 

In  subtracting  the  proper  background  from  the  raw  Auger  data 
one  normally  utilizes  low  order  polynomials  for  weakly  varying 
backgrounds  like  that  shown  in  Fig.  la.  For  our  present  purposes 
only  a  slight  linear  background  correction  was  necessary  and  the 
result  of  such  a  background  subtraction  is  shown  in  Fig.  1b.  The 
criteria  for  a  properly  corrected  background  are:  (1)  that  the 
region  of  the  corrected  spectrum  with  energies  above  the  Auger 
threshold  ( >280  eV  in  Fig.  1b)  should  be  flat  and  at  the  zero 
basbline  and  (2)  that  the  featureless  region  below  the  structure  in 
the  background-corrected  data  (<230  eV  in  Fig.  1b)  should 
accurately  match  in  shape  the  corresponding  region  in  the  loss 
function.  The  shape  similarity  In  this  region  is  evident  from  the 
comparison  .shown  in  Fig3.  1c  and  Id.  If  these  criteria  are  not  met 
then  one ’of  the  basic  assumptions  concerning  the  deconvolution 


procedure  Is  not  fulfilled  and  a  proper  deconvolution  cannot  be 
made . 

Since  the  true  loss  function  is  noU  available ,  one  other 
assumption  Is  necessary  for  proper  loss  correction  by 
deconvolution.  We  must  assume  that  the  general  shape  of  the  model 
function  is  the  same  as  the  true  loss  function,  the  only  possible 
difference  being  the  relative  intensities  of  the  loss  features  with 
respect  to  that  of  the  elastic  peak.  In  the  commonly  used 
procedure  (31),  the  losses  are  stripped  away  from  the  elastic  peak, 
the  elastic  peak  is  normalized  to  a  unit  area  and  then  the  losses 
are  replaced  scaled  by  an  adjustable  parameter.  The  deconvolution 
is  then  performed  interactively  adjusting  the  scaling  parameter  to 
achieve  a  flat  and  zero  baseline  behavior  in  the  low-energy  region 
beyond  the  characteristic  features. 

The  result  of  this  procedure  for  the  data  of  Figs.  1b  and  1e 
is  shown  in  Fig.  2.  Careful  attention  to  the  criteria  just 
described  results  in  corrected  spectra  which  are  very  reproducible, 
probably  to  within  _*  20t,  In  the  low-energy  region  of  the  Auger 
line  where  the  deconvolution  procedure  Is  most  susceptible  to 
error.  Any  distortions  between  the  results  given  in  Fig.  2  and  the 
"true"  Auger  spectrum  are  minimal  in  the  threshold  region  and 
progressively  Increase  through  the  low-energy  tail. 

Significant  angular-dependent  lineshape  variations  have  been 
found  in  the  main  body  of  the  Auger  spectrum  for  HOPG  and  SCG 
samples  (39).  In  order  to  better  approximate  the  graphite  spectrum 
in  terms  of  one-electron  models  which  are  angle  i ntegi For 
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nature,  we  have  chosen  to  use  the  POCO  graphite  results  of  Fig.  1 
in  our  subsequent  analyses. 


The  Fermi  level  noted  by  the  FL  lin»  in  Fig.  2,  was  accurately 
established  in  a  separate  experiment  in  which  he  X-ray  excited 
Auger,  the  core  level,  and  valence  spectra  were  obtained 
simultaneously.  The  Fermi  level  expressed  on  a  KE  scale  is  simply 
the  measured  C(1s)  binding  energy  of  28H.6  eV.  An  example  of  this 
data  is  shown  in  Fig.  3.  These  results  permit  an  accurate 
calibration  of  the  Auger  kinetic  energy  with  respect  to  the  Fermi 
energy. 

It  should  be  noted  with  respect  to  the  corrected  Auger 
lineshape  of  Fig.  2  that  similarly  obtained  results  reported 
earlier  by  Smith  and  Levenson  {23}  show  a  considerably  narrower 
line  with  no  apparent  low-energy  structure.  However,  the 
background-corrected  raw  data  and  elastic/near-elastic  backscatter 
spectra  shown  in  the  Smith  and  Levenson  paper  (Figs.  1  and  2  of 
ref.  23)  do  not  satisfy  the  shape  criterion  in  the  low-energy  tail 
region  as  we  have  just  discussed.  The  result  of  this  mismatch  is 
that  the  deconvolution  forces  the  corrected  spectrum  to  go  negative 
in  the  low-energy  region.  Using  our  spectra  in  Figs.  1b  and  1e, 
the  deconvoluted  results  of  Smith  and  Levenson  can  only  be 
approximately  obtained  by  a  serious  distortion  of  the  deconvolution 
procedure.  Apparently,  Smith  and  Levenson  truncated  their 
corrected  spectrum  at  the  point  at  which  it  went  negative  (due  to 
improper  removal  of  the  featureless  tail)  giving  the  impression  of 
a  narrower  line  with  no  intensity  in  the  low-energy  region. 
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Because  of  the  possibility  of  distorted  lineshapes  similar  to 
that  experienced  by  Smith  and  Levenson  resulting  from  the  data 
reduction  procedures,  we  have  taken  particular  care  to  evaluate  the 
reliability  of  the  experimental  data  and  the  data  manipulation 
procedures  presented  here.  We  have  obtained  corrected  Auger 
llneshape  data  independently  from  three  different  graphite  samples, 
three  different  spectrometers  and  at  two  separate  laboratories. 

This  procedure  resulted  in  very  similar  lineshapes  for  similarly 
prepared  POCO  graphite  samples.  The  relative  intensity  of  the 
feature  in  the  low-energy  region  at  about  245  eV  (which  is  the  most 
sensitive  to  the  data  manipulation  techniques  for  all  samples)  was 
found  to  vary  in  our  experiments  by  about  20J,  but  the  feature  was 
always  present.  With  this  relative  accuracy  in  mind,  in  the 
remainder  of  this  paper  we  will  attempt  to  clarify  the  origins  of 
the  graphite  llneshape. 

IV.  DISCUSSION 

In  attempting  to  characterize  the  Auger  llneshape  for  graphite 
in  terms  of  its  known  electronic  properties  it  is  helpful  to  begin 
as  simply  as  possible  and  allow  the  disagreement  between  the  model 
characterization  and  the  experimental  result  to  guide  further 
sophistication.  For  this  approach,  the  simplest  model  involves  the 
one-ielectron  approximation  with  the  inclusion  of  the  Auger  matrix 
elements  based  on  similar  transitions  in  rare-gas  atoms.  We  begin 
the  present  discussion  with  such  a  comparison  and  continue  by  using 
the  nature  of  the  disagreement  as  a  guide  to  the  inclusion  of 
"many-electron"  effects  such  as  Initial-state  screening,  shake 


* 


-  11  -  <1 

phenomena,  hole-hole  correlation,  and  dynamic  final-state 
screening. 

A.  One-electron  Approach 

Murday,  et  al .  {21}  have  deduced  from  experimental  data  the 
graphite  partial  density  of  states  (PDOS)  components:  Up,  op  and 

o  •  The  p  electrons  contribute  to  both  the  »  and  o  bands.  The  o„ 

PDOS  was  obtained  from  the  XPS  valence  band  (38)  where  the 
intensity  is  primarily  determined  by  s  symmetry  (i.e.,  XPS 
intensity  is  os  ♦  1/32(op  ♦  ir  )){21).  The  op  and  x  PDOS  were 

obtained  from  angular-dependent  XES  data  {40}:  the  op  PDOS  was 

obtained  from  data  at  80°  take-off  angle  and  the  irp  from  5°  data. 

The  relative  areas  of  the  individual  PDOS  were  normalized  to 
2 

provide  an  o3opit  configuration.  This  total  empirical  DOS  has  been 

found  to  agree  reasonably  well  with  theoretical  calculations  {41- 
43)  except  for  the  overall  bandwidth.  Murday  et  al .  {21},  forced 
the  bandwldths  to  agree  by  arbitrarily  truncating  the  PDOS. 

Such  truncation  is  reasonable  since  the  empirical  procedures 
for  determining  PDOS  all  involve  experimental  results  from 
excitation  or  deexcitation  techniques  that  are  not  in  themselves 
one  electron  in  nature.  These  techniques  invariably  show  bands 
that  are  broader  than  expected  theoretically,  presumably  due  to  the 
effects  of  intrinsic  processes,  and  truncation  has  become  a  common 
procedure  to  correct  for  such  effects  {21,38,44}.  However,  since 
we  are  trying  to  model  a  many-electron  process  by  manipulating  the 
results  of  other  many-electron  processes,  it  seems  more  reasonable 
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to  utilize  the  data  itself  rather  than  perform  an  arbitrary 
truncation.  Therefore,  we  show  the  results  of  the  Murday,  et  al . 
procedure  without  data  truncation  as  th'e  solid  lines  in  Fig.  4. 

2 

Figure  4a  contains  the  og  and  Op  components,  scaled  for  a  ogOp  n 

electron  configuration,  of  the  PDOS  as  obtained  from  the  XES  and 
XPS  data;  Fig.  4b  contains  the  o  and  it  PDOS  (o-o  +o  );  and  Fig.  4c 

3  p 

is  the  total  DOS  (o  ♦  it)  .  The  dotted  section  of  Fig.  4c 
illustrates  the  truncated  results  published  by  Murday  et  al .  (21). 
The  XES  and  XPS  spectra  leading  to  the  PDOS  in  Fig.  4  were 
originally  referenced  to  the  graphite  Fermi  energy  and  the  energy 
scale  shown  in  Fig.  4  was  obtained  by  adding  the  graphite  C(1s) 
binding  energy  (284.6  eV)  to  this  Fermi  energy.  Thus,  both  the 
empirical  PDOS  data  and  the  experimental  Auger  lineshape  of  Fig.  2 
are  placed  on  a  common  absolute  energy  scale  and  the  subsequent 
modeling  will  be  done  without  altering  these  energy  scales. 

As  mentioned  earlier,  the  model  Auger  lineshape  A(E)  in  the 
one-electron  approximation  can  be  generated  by  self-folding  the 
DOS.  In  terms  of  the  three  PDOS  components  shown  in  Fig.  4a,  this 
procedure  results  in  the  expression, 

A(E)  ■  F>kSS°S*0S  +  2Pksp(0S*0p  +  09*1rp)+Pkpp(0p*0p+20p,,irp+1Ip#1,p)  (2) 

where  the  o3*og,  etc.  terms  indicate  a  fold  of  the  PDOS  components, 

o’o.  -  /o„(E)o  (E-e)de,  (3) 

and  the  Pkn  factors  are  atomic  Auger  matrix  elements  normalized 

per  electron.  These  atomic  matrix  elements  are  assumed  to  be  the 
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same  as  those  determined  from  both  experiment  and  theory  for  atomic 


neon  {45} : 


A  /A 
kss  kpp 


0.0S9,  Aks|>  /»kpp  -  031. 


A  plot  of  the  relative  contributions,  A^^,,  as  a  function  of  Z, 

t 

the  nuclear  charge,  reveals  that  the  ratios  (11  -  ss  and  3p),  are 

remarkably  constant  from  atom  to  atom  {8}.  From  the  expressions, 

Akss  -  ^kss  and  Aksp  -  ^Pksp 
Akpp  ^Pkpp  Akpp  ^6Pkpp 

the  relative  atomic  Auger  matrix  elements  per  electron,  P^gg'^sp 
: P.  -  8:0. 5:1.0,  can  be  obtained.  Although  these  same  matrix 

kpp 

elements  were  utilized  by  Murday,  et  al.  {21},  the  factor  of  two  in 

front  of  the  o  *n  term  in  eq.  2  was  Inadvertently  omitted  in  their 
P  P 

work.  Thus,  the  o  *rr  Intensity  in  their  Fig.  2  is  a  factor  of  two 
P  P 

too  small. 

We  can  now  use  eq.  2  and  the  nontruncated  PDOS  functions  of 
Fig.  4a  to  generate  the  one-electron  model  function  for  the  Auger 
llneshape  which  is  shown  as  the  dotted  curve  in  Fig.  5.  The  dashed 
portion  shows  the  effect  of  truncating  the  PDOS  {21}.  The 
experimental  result  from  Fig.  2  is  shown  a3  the  solid  curve  and  it 
is  apparent  from  this  comparison  that  the  one-electron  model  of  the 
Auger  llneshape  significantly  differs  from  the  experimental  result 
in  several  respects  regardless  of  whether  a  truncated  or 
nontruncated  PDOS  is  used.  The  model  llneshape  'a  appreciably 
narrower  than  the  experimental  result  with  intensity  missing  in 
both  the  threshold  region  n^ar  280  eV  and  in  the  area  below  the 


principal  maximum  at  about  260  eV.  Thus,  it  is  clear  that 
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processes  other  than  those  incompassed  by  the  one-electron  model 
play  a  significant  role  in  determining  the  graphite  Auger 
lineshape . 

B.  Many-electron  Effects 

The  one-electron  approximation  assumes  that  the  Auger 
transition  begins  with  an  excited  core  state  (core  hole)  and  ends 
with  two  holes  in  valence  levels  along  with  an  Auger  electron  in 
the  continuum.  During  this  entire  process,  the  valence  electronic 
levels  are  considered  to  be  frozen  in  their  ground-state 
configurations,  which  would  imply  that  the  electrons  of  the  system 
are  noninteracting.  Of  course,  if  the  letter  true  there  would  be 
no  Auger  transition  so  it  is  obvious  that  the  approximation  is 
going  to  break  down  and  the  only  question  is,  "how  does  this  break 
down  manifest  itself  in  the  Auger  lineshape  and  to  what  relative 
magnitude?"  We  will  structure  our  discussion  of  many-electron 
effects  as  if  they  were  separable  in  their  contribution  to  the 
Auger  lineshape  which  is  somewhat  artificial  since  the  effects  are 
all  part  of  the  same  phenomenon.  The  various  effects  will  be 
discussed  ac'ording  to  their  contributions  to  the  different  regions 
of  the  -pectrum  including  those  which  principally  affect  the 
intensity  of  the  self-fold  components,  those  which  give  rise  to 
intensity  above  threshold,  and  those  which  contribute  to  low-energy 
intensity. 

1.  Static  Initial-state  Screening 

The  static  screening  response  of  the  valence  electrons  to  the 
presence  of  the  core  hole  in  the  Auger  initial  state  {30, 46}  has 
been  shown  to  contribute  to  the  intensity  of  the  various  self-fold 

t 


15 


<• 

components.  The  term  "static"  differentiates  this  aspect  of  the 
screening  response  from  the  dynamic  shake  phenomena  that  will  be 
discussed  below.  This  distinction  has  been  discussed  recently  by 
Sawatzky  {47}  with  respect  to  core-hole  creation  in  photoelectron 
spectroscopy.  If  one  slowly,  or  adiabatically ,  creates  a  core  hole 
in  a  solid-state  atom  the  electronic  states  will  adjust  themselves 
in  an  attempt  to  lower  the  energy  of  the  core-hole  state.  The 
static  screening  response  amounts  to  the  situation  where  the  core¬ 
hole  creation  results  in  the  fully  relaxed  state  with  the 
electronic  states  polarized  as  if  they  were  responding  to  the 
presence  of  the  Z  +  1  impurity  3tate. 

Static  screening  has  been  discussed  for  Auger  transitions  by 
Ramaker  and  co-workers  {26,46,48},  and  by  Jennison  {30}  and  can  be 
illustrated  through  the  use  of  the  final-state  rule  (FS  rule).  The 
FS  rule  indicates  that  the  total  intensities  of  the  various 
symmetry  components  making  up  the  Auger  lineshape  (ss,  sp,  etc.) 
are  determined  by  their  electronic  configuration  appropriate  to  the 
statically  screened  initial-state  core  hole,  while  the  shape  of 
each  contribution  is  determined  by  the  appropriate  final-state  DOS 
{46}.  If  we  further  assume  that  the  final-state  holes  are 
noninteracting,  then  according  to  the  FS  rule  the  ground  state  DOS 
can  ‘be  substituted  for  the  final  DOS  as  was  done  by  Murday,  et  al . 
{21}. 

The  statically  screened,  initial-state  configuration  has  been 
determined  recently  from  a  parameterized  109  atom  cluster  LCA0-M0 
calculation  by  Dunlap,  et  al .  {26}  under  the  assumption  that  only 
the  unfilled  »r  band  makes  a  significant  contribution.  It  was  found 
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that  the  local  electron  density  In  the  it  band  increased  by  0.55 

2  1  55 

electrons  to  a  o  o  II  configuration.  However,  recent 
s  p 

calculations  by  Binkley  {49}  on  benzene  and  pyridine  suggest  that 
the  screening  response  of  the  other  components  of  the  valence  band 
are  also  important.  Benzene  with  a  core  hole  at  one  carbon  site 
and  the  ground  state  electron  density  at  the  nitrogen  site  in 
pyridine  give  essentially  the  same  result  indicating  that  the  it 
screening  charge  is  increased  by  about  0.5*1  (in  agreement  with 
Dunlap,  et  al .  {26])  and  that  the  o  and  o  states  pick  up  about 

P  3 

0.23  and  0.52  electrons,  respectively  as  a  result  of  static, 

1 nl tl al -state  screening. 

These  variations  in  the  i ni ti al -state  electron  occupation 

result  in  only  minor  changes  in  the  model  Auger  lineshape  as 

illustrated  in  Fig.  6.  Here,  the  solid  curve  includes  only  the  ir- 

level  screening  of  Dunlap,  et  al.  {26)  and  the  dotted  curve 

Includes  screening  from  all  the  band3  using  the  Binkley  {*19} 

parameters.  The  level  of  agreement  with  the  experimental  spectrum 

is  only  slightly  changed  in  either  case. 

The  model  lineshapes  in  Fig.  6  were  obtained  utilizing  the  FS 

rule  as  stated  above,  which  is  known  to  break  down  within  a  few  eV 

of  the  Fermi  level  {21 }  suggesting  that  the  difference  near  the  top 
* 

of  the  lineshape  between  the  model  and  experiment  might  result  from 
this  breakdown.  However,  our  model  lineshape  was  obtained  by 
folding  expirically  determined  PDOS  from  the  XES  spectra  and  if  the 
FS  rule  is  seriously  breaking  down  for  the  Auger  lineshape,  it 
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would  similarly  break  down  for  the  XES  lineshape.  Thus,  initial- 
state/final-state  nonorthogonality  effects  near  the  Fermi  level 
have  been  Included  implicitly  in  the  moc?el.  In  addition,  the 
calculations  by  Dunlap,  et  al .  {26}  indicate  that  these 
nonorthogonality  effects  are  relatively  small  for  graphite. 

2.  Dynamic  Ini ti al -state  Screening 

As  we  have  shown,  static,  initial-state  screening  cannot 
account  for  the  differences  between  the  model  lineshape  and 
experiment  for  the  region  of  the  Auger  spectrum  Ju3t  below  the 
Fermi  level.  However,  there  are  several  dynamic  screening  effects 
involving  the  initial  state  which  can  give  rise  to  structure  in 
this  region.  For  example,  the  sudden  creation  of  a  core  hole  can 
lead  to  shakeup  processes  which  leave  the  system  in  a  more 
energetic  initial  state  than  a  statically  screened  core  hole.  If 
this  excited  situation  remains  local  to  the  core-hole  site  for 
sufficient  time,  the  Auger  transition  can  utilize  electrons  n  the 
excited  state  and  produce  Intensity  at  higher  energies  than 
expected  on  the  basis  of  a  static  screening  model. 

Shakeup  structure  of  this  kind  is  well  known  in  gas-phase 
molecules  {50},  transition  metals  {51},  semiconductors  and 
insulators  {52}.  Such  excited  states  can  readily  be  produced  by 
electron  and  off-resonant  photon  excitation  near  the  core-hole 
excitation  threshold  and  specific  excited  states  can  be  produced  by 
resonant  photon  excitation.  Graphite  is  known  to  have  a  prominent 
resonant  excitation  structure  just  above  the  Fermi  level  resulting 
from  a  core-exci toni c  state,  i.e.,  an  excited  state  consisting  of  a 
conduction-band  electron  bound  to  the  excited  core  hole.  The 
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existence  of  this  state  was  established  by  Mele  and  Rltsko  {53) 
using  electron  energy  loss  spectroscopy  with  80  keV  Incident 
electrons.  In  their  spectra  the  core  excltonic  state  is  seen  as  a 
sharp  level  centered  at  about  1  eV  above  the  Fermi  level  with  a 
full  width  at  half  maximum  of  about  1.0  eV. 

If  this  core-excitonic  state  were  Involved  in  the  subsequent 
Auger  transition,  then  we  would  expect  to  see  intensity  in  the 
high-energy  region  of  the  Auger  spectrum  in  electron-excited 
results  but  not  in  off-resonant ,  photon-excited  measurements.  In 
fact,  we  find  identical  lineshapes  using  electron  excitation  and 
excitation  by  Mg(Ka)  photons  (1253.6  eV).  Thus,  it  is  apparent 
that  direct  excitation  into  the  core-excitoni c  state  does  not 
result  in  appreciable  Auger  intensity. 

The  reason  that  resonant  excitation  into  the  core-excitonic 
3tate  does  not  contribute  significant  Auger  intensity  is  no  doubt 
due  to  the  short  lived  existence  of  an  electron  in  this  state.  From 
the  lifetime  broadening  of  the  C  ( 1 s )  core  state  {5*0  and  the  core- 
excitonic  state  153)(0.06  eV  and  1.0  eV,  respectively),  we  estimate 
that  the  lifetime  of  the  former  is  about  17  times  the  latter. 

Thus,  there  is  only  a  small  probability  that  this  directly 
populated  core-excitonic  state  would  be  occupied  during  the  Auger 
transition . 

Van  Attekum  and  Wertheim  {55}  have  shown  that  the  C(1s)  XPS 
spectral  line  shape  of  graphite  is  distorted  due  to  the  shakeup  of 
a  valence  electron  into  an  excl toni c-1 ike  level  near  the  Fermi 
level  during  the  core  excitation.  This  process  differs  from  the 
direct  excitation  in  that  the  resultant  Auger  initial  state,  which 
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we  will  term  a  valence/core  exciton,  would  contain  two  positive 
holes,  one  in  the  core  level  and  one  in  a  valence  level,  with  one 
electron  in  the  excitonlc  level.  (The  direct  excitation  process 
has  a  missing  core  electron  and  an  electron  occupying  the  excitonlc 
level).  Their  analysis  suggested  a  lifetime  broadened  width  ior 
this  excited  condition  of  approximately  0.1  eV  indicating  a 
lifetime  an  order  of  magnitude  longer  than  the  directly  populated 
core  exciton  (1.0  eV)  and  one  approaching  that  of  the  core  hole 
(0.06  eV).  It  is  probable  that  this  valence/core  exciton  doe3  make 
a  contribution  to  the  threshold  region  of  the  graphite  Auger 
spectrum.  Moreover,  since  the  fluorescent  yield  for  graphite  154) 
Indicates  a  photon  decay  lifetime  width  of  only  0.0002  eV,  the 
valence/core  exciton  is  not  expected  to  contribute  to  the  XES 
spectrum . 

Inclusion  of  the  valence/core  exciton  results  in  significant 
lineshape  changes  since  it  gives  rise  to  intensity  in  the  threshold 
region  just  below  the  Fermi  level  where  intensity  is  missing  in  the 
model  spectra  in  Figs.  5  and  6.  This  can  be  seen  in  Fig.  7  where 
we  have  included  the  electron  density  arising  from  shakeup  into  the 
valence/core  excitonlc  state  as  a  delta  function  at  the  Fermi  level 
and  varied  its  electron  occupancy  to  obtain  a  "best  fit"  with  the 
leading  edge  of  the  experimental  result.  We  have  assumed  that  the 
excited  electron  in  the  valence/core  excitonlc  state  has  p  symmetry 

and  have  used  the  appropriate  Auger  matrix  elements  discussed 

> 

earlier.  The  value  for  the  effective  electron  occupancy  in  the 
valence/core  excitonlc  state  obtained  by  this  procedure  is  0.27 
electrons  which  appears  reasonable  compared  with  the  intensity  in 


the  lower  shoulder  of  the  core  XPS  peak  {55}.  A  full  account  of 
the  origin  and  modeling  of  this  feature  is  given  elsewhere  {56}. 

The  additional  intensity  in  the  high-energy  region  of  the 
spectrum  shown  in  Fig.  7  takes  two  forms;  (1)  a  sharp  feature 
located  near  the  Fermi  level  resulting  from  the  Auger  transition 
which  leaves  two  holes  in  the  valence/core  excitonic  state  and  (2) 
a  broader  structure  located  below  the  Fermi  energy  resulting  from 
the  transition  which  places  one  hole  in  the  valence/core  exci*'  nic 
level  and  one  in  the  valence  band.  The  fact  that  the  experimental 
Auger  spectrum  shows  no  sharp  feature  near  the  Fermi  energy  simply 
means  that  the  effective  occupancy  of  the  valence/core  excitonic 
state  is  not  sufficient  to  give  rise  to  a  significant  probability 
of  it  being  doubly  occupied.  Simple  statistical  arguments  suggest 
that  an  average  effective  occupancy  of  0.27  electrons  in  this  state 
would  lead  to  an  intensity  ratio  of  about  5?  between  the  Auger 
feature  resulting  from  both  holes  in  the  valence/core  excitonic 
state  compared  to  one  hole  in  this  state.  Thus,  it  is  not 
surprising  that  this  feature  is  not  seen  above  the  noise  level  in 
the  experimental  spectrum. 

It  is  interesting  to  note  that  in  intercalated  graphite  with  a 
donor  lntercalant,  e.g.,  alkali  metals,  the  Fermi  level  position 
inepeases  as  some  of  the  normally  empty  conduction-band  states  are 
filled  {26}.  This  process  gives  rise  to  an  increased  occupancy  of 
the  core  excitonic  Auger  initial-state  and  a  sharp  feature  of 
significant  intensity  near  the  Auger  threshold.  Furthermore,  the 
lifetime  of  this  core-excl toni c  screening  charge  is  long  since  it 
i 3  part  of  the  static  screening  process  and  lies  below  the  Fermi 
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level.  The  Intensity  of  this  cor e-e xci toni c  contribution  is  found 
to  increase  as  the  square  of  the  intercalant  concentration ,  whereas 
the  intensity  of  the  structure  which  would  correspond  to  that 
enhanced  in  Fig.  7  only  increases  in  a  linear  fashion  {27}.  This 
is  to  be  expected  on  the  basis  of  the  probability  of  double-hole 
occupancy  of  the  core  excitonic  state  relative  to  single-hole 
occupancy . 

From  the  comparison  of  Fig.  7,  we  can  see  that  dramatic 
Improvement  is  obtained  in  the  near  Fermi-level  portion  of  the 
model  by  including  the  effect  of  the  dynamic  1 ni t i al -state  Auger 
processes.  However,  significant  discrepancies  remain  in  the  area 
below  the  principal  peak  and  it  is  apparent  that  other  many- 
electron  effects  must  be  taken  into  account. 

3.  Final  State  Hole-hole  Interaction 

The  Auger  transition,  being  two-electron  in  nature,  causes  a 
unique  man y-el ect ron  effect  which  significantly  distorts  the 
lineshape  compared  to  that  expected  from  a  one-electron  model.  The 
effects  of  hole-hole  Interaction  on  Auger  llneshapes  has  been 
modeled  for  simple  systems  by  Cini  {35,36}  and  Sawatzky  {37}.  The 
holes  tend  to  remain  spatially  localized  after  the  Auger  transition 
giving  rise  to  a  net  Coulomb  repulsion  characterized  by  the 
effective  energy  value  Ueff.  Cini  has  given  the  following 

expression  {35,36}  for  the  distortion  caused  by  localization  of  the 


Auger  final-3tate  holes 
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where  N(E)*N(E),  the  self-fold  of  the  one-electron  DOS,  N(  E) , 
represents  the  undistorted  Auger  lineshape  in  the  one-electron 
approximation  of  eq.  2,  and  I^( E)  is  the  Hilbert  transform  of 

N(E)*N(E).  The  Cini  expression  was  derived  from  a  many-electron 
calculation  on  a  single  filled  band. 

As  was  demonstrated  by  Cini  {36},  the  effect  of  the  hole-hole 
interaction  is  to  shift  spectral  intensity  toward  the  low-energy 
end  of  the  lineshape  as  the  value  of  Ueff.  is  increased.  The  energy 

positions  of  the  top  and  bottom  of  the  spectrum  do  not  change  from 
their  one-electron  values.  As  the  value  of  becomes  greater 

than  the  width  of  the  original  valence  band,  a  discrete  state  can 

be  split  off  below  the  main  Auger  structure.  This  discrete 

structure  is  expected  to  be  much  narrower  and  more  intense  than 

that  of  the  main  portion  of  the  Auger  spectrum. 

In  attempting  to  model  the  effect  of  the  hole-hole  interaction 

in  grapmte  using  the  Cini  expression,  we  must  first  address  the 

fact  that  we  are  not  dealing  with  a  single,  filled-band  system.  In 

fact,  in  graphite  we  have  two  bands  as  shown  in  Fig.  Mb  (the  0  and 

x  bands)  both  of  which  are  half  filled.  In  addition,  the  atomic 

orbital  composition  of  the  0  band  varies  with  position  in  the  band 

as  shown  in  Fig.  Ma.  To  our  knowledge,  no  theoretical  expression 
* 

is  presently  available  corresponding  to  eq.  5  for  half-filled 
bands,  although  Cini  and  others  have  discussed  the  case  of  a  small 
number  of  empty  states  in  an  otherwise  filled  band  {57,58}. 
Significantly  unfilled  bands  complicate  the  situation  by 
Introducing  an  increased  screening  response  in  such  systems.  In 
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the  absence  of  an  adequate  theoretical  treatment  for  this  case,  we 
will  assume  that  the  mixing  of  the  empty  states  in  the  o  and  u 
bands  is  small  and  will  only  affect  the  value  of  Uef.f  in  eq.  5. 

This  assumption  seems  reasonable  since  the  antibonding  portion  of 
the  o  band  is  split  off  from  the  bonding  portion  and  we  have  shown 
previously  that  as  long  as  is  small  relative  to  this 

separation,  the  bonding  and  antibonding  bands  do  not  mix  (59). 

This  assumption  may  not  be  valid  for  the  ir  band,  but  the  excellent 
agreement  between  the  one-electron  model  with  dynamic  initial-state 
screening  effects  and  the  experimental  Auger  lineshape  near  the 
Fermi  level  {56},  indicates  that  localization  effects  here  are 
neglible  anyway. 

A  problem  with  using  the  Cini  expression  on  a  filled 
degenerate  two  band  system  arises  from  the  possibility  that  the 
bands  could  mix  giving  rise  to  cross  terms  in  the  distortion 
expression.  Cini  has  shown  formally  that  the  extension  of  his 
theory  to  degenerate  orbitals  and  bands  can  be  solved  exactly,  but 
the  calculations  and  results  are  very  complicated  {36}.  He  has 
suggested  that  if  the  solid  does  not  distort  the  spherical  symmetry 
of  the  atom  significantly,  such  as  in  metals,  the  equations  for  the 

23  4  1 

different  total  angular  momenta,  decouple  so  that  each  L 
raulttplet  component  can  be  treated  as  an  independent  band  using  eq. 
5  with  a  different  U  parameter  for  each.  Previously,  Weightman 

and  Andrews  {60,61}  analyzed  their  Auger  spectra  of  transition 
metals  and  alloys  with  this  approximation  and  obtained  good 
agreement  with  experiment. 
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The  spherical  symmetry  appproximation  is  not  valid  for 
graphite,  or  for  that  matter,  for  any  covalently  bonded  solid,  but 
we  can  make  a  similar  approximation  utilizing  the  local  D^,  planar 

symmetry  around  each  carbon  atom.  Utilizing  this  approximation 

2s  +  i 

each  multiplet  decouples  in  a  manner  similar  to  that  used  by 

Ramaker  {62)  for  a  molecular  orbital  cluster  calculation  for  Si02. 

However,  the  number  of  different  multiplets  and  U  parameters 

which  arise  becomes  large  and  the  essential  physics  is  difficult  to 
ascertai n . 

We  prefer  a  more  intuitive  approach  which  we  believe  is  valid, 
since  the  magnitude  of  the  correlation  effects  in  graphite  are 
certainly  much  less  than  in  either  S102  or  in  the  transition  metal 

alloys  {60,61).  We  assume  that  multiplets  coming  from  the 
individual  o*o,  o*u  and  it*ir  folds  have  a  common  Uef{,  value,  thus, 

reducing  the  number  of  parameters  to  three.  Furthermore  we  assume 
that  the  relative  separation  in  energy  between  the  o*o,  o*m  and  tt*h 
contributions  is  large  compared  to  so  that  the  three 

contributions  can  be  considered  separately.  The  effect  of  the 
Auger  matrix  elements  will  be  handled  by  projecting  the  symmetry 
components  onto  the  Cini  distortion  function  and  multiplying  the 
result  by  the  appropriate  matrix-element  values.  As  an  example, 
for  the  o  band  we  first  write  the  Cini  distortion  function  from  eq. 
5  a3 
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where  N  *N  is  the  total  o  band  self-fold 
o  o 
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Therefore,  the  three  Cini  distorted  Auger  contributions  are  given 
by  the  expressions 
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terms  are  the  relative  matrix 


elements  given  earlier.  We  now  have  just  three  components  with 
three  different  values  of  and  the  model  lineshape  is  given  by 


A(  £)  -  A  (  E)  ♦  A  (E)  ♦  A  ( E)  .  (11) 

oo  on  IT  IT 

Using  this  procedure,  we  are  able  to  produce  the  "best  fit" 

data  shown  a3  the  dashed  curve  in  Fig.  8.  The  dotted  curves  A,  B 

and  C  are,  respectively,  the  contributions  A  ,  A  and  A  Trom 

oo  ott  uir 

eqs.  8-10,  and  D  is  the  valence/core  exciton  contribution  (56). 

This  particular  model  Auger  lineshape  was  obtained  with  Uef.j.  values 

of  2.2,  1.5,  and  0.6  eV  for  the  o*o,  o*n,  and  n*n  self-folds, 
respectively.  The  near  zero  value  of  Ugj-j-  for  the  ir*n  self-fold  is 

consistent  with  previous  data  for  benzene  {10}  where  Ue(lf.  for  the  n 


transitions  reflects  delocalization  to  the  dimensions  of  the 
molecule.  We  note  that  there  has  been  considerable  improvement 
relative  to  the  comparison  of  Fig.  5  in  the  threshold  region  and  in 
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the  area  of  the  principal  maximum.  However,  significant 
discrepancies  remain  in  the  region  below  245  eV. 

The  general  criterion  for  atomic  localization  obtained  from 
the  Cini-Sawatzky  approach  for  single  metallic  bands  depends  on  the 
relative  magnitude  of  the  bandwidth  W  compared  to  the  value  of 
Ueff*  Here  w  refers  to  the  total  bandwidth  (bonding  plus 

antibonding)  and  is  a  measure  of  the  strength  of  the  covalent 

interaction.  For  U  <<  W  the  final  states  are  delocalized  while 

ef  f 

for  u  >>  W  the  final  states  are  strongly  localized.  Similar 

conclusions  have  been  reached  for  covalent  systems  in  terms  of 
atomic,  bonding  and  group  orbital  contributions  by  Dunlap,  et  al . 
(63). 

The  empirical  values  of  Ueff.  obtained  for  the  self-folds  from 

the  model  fit  in  Fig.  8  appear  to  be  reasonable  based  on  a  recent 

application  ot  these  criteria  to  graphite  (69).  The  screening  in 

graphite  would  certainly  be  expected  to  result  largely  from  n-n* 

mixing  and  this  screening  would  be  most  effective  for  ir-band  holes. 

As  mentioned,  structure  in  the  low-energy  region  of  the 

lineshape,  where  no  intensity  is  expected  on  the  basis  of  the  one- 

electron  model,  could  presumably  be  obtained  by  simply  Increasing 

the'  value  of  U  to  the  extent  of  creating  a  discrete  state  below 
erf 

the  main  portion  of  the  line.  However,  estimates  of  from 

LCAO-ZDO  calculations  using  various  models  for  the  effect  of 
screening  on  U e f r ,  suggest  that  values  for  Uef  f  are  not  sufficient 

to  give  rise  to  strong  localization.  Thus,  the  broad  feature  in 
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Fig.  8  centered  at  about  2^0  eV  must  be  due  to  effects  other  than 
J.ole-hole  localization. 

*».  Dynamic  Final-state  Screening 

In  the  present  section  we  briefly  outline  the  effects  on  the 
Auger  lineshape  of  the  dynamic  aspects  of  final-state  screening. 
Auger  transitions  which  themselves  leave  the  system  in  an  excited 
state  will  produce  structure  at  energies  below  that  expected  from 
the  one-electron  picture.  These  result  from  the  so  called  "shake" 
phenomena  in  the  final  state.  Cini  and  D'Andrea  { 64 }  have  recently 
discussed  the  effects  of  dynamic  screening  of  the- hole-hole 
Interaction  which  can  also  lead  to  structure  below  that  expected 
from  the  one- electron  model. 

The  structure  near  2M0  eV  may  indicate  that  dynamic  final- 
state  screening  is  contributing  appreciably  to  the  overall 
lineshape.  The  Auger  transition  proceeds  from  a  screened  core-hole 
3tate  to  that  of  a  neutral  core  and  a  two-hole  state  in  the  valence 
levels.  From  the  standpoint  of  the  core  state,  thi3  transition  is 
just  the  reverse  of  the  core  excitation  process  and  should  give 
rise  to  the  same  array  of  dynamic  screening  features  seen  in  the 
initial  state. 

The  screening  associated  with  the  creation  of  a  core  hole  in 
graphite  has  been  extensively  studied  by  a  number  or  experimental 
techniques.  Core-level  XPS  spectra  show  characteristic  loss 
structure  associated  with  pl^sraon  excitation  (38}  very  similar  to 
that  shown  in  the  backscatter  spectrum  of  Fig.  1e.  These 
techniques,  however,  yield  information  regarding  the  extrinsic  loss 
processes  as  well  as  the  intrinsic  dynamic  screening  processes  and 
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it  is  difficult  to  delineate  the  contribution  of  the  latter  with 
reliable  accuracy  (65,66}. 

Clear  evidence  for  the  intrinsic  processes  in  graphite  has 
been  obtained'  from  XES  where  there  are  no  transport  electrons  to 
cause  extrinsic  loss  processes  (67).  Furthermore,  core-level  XPS 

for  gas-phase  benzene,  where  extrinsic  loss  processes  should  be 

•  ) 

greatly  suppressed,  shows  loss  structure  very  similar  to  that  of 
graphite  (68-70). 

The  dynamic  aspects  of  the  response  to  the  creation  of  the 
initial-state  core  hole  consist  of  a  narrow  excitation  at  a  loss 
energy  of  about  7  eV  (associated  with  the  creation  of  a  plasmon 
Involving  the  ir-ir*  states  {71}),  a  broad  excitation  centered  at  a 
loss  energy  of  approximately  27  eV  (attributed  to  a  plasmon 
creation  Involving  a  combination  of  o,  and/or  it,  and  o*  states 
(71)),  a  broad  featureless  tail  related  to  interband  transitions, 
and  a  distortion  of  the  symmetry  of  the  core-level  photoemission 
line  associated  with  shake-up  into  the  core  excitonlc  state  (55). 
With  the  exception  of  the  core-excitonlc  state,  all  of  these 
characteristics  of  the  core  hole  screening  processes  (Including 
conduction/valence-band  electron/hole  creation)  can  be  seen  in  the 
backscatter  loss  data  for  graphite  shown  in  Fig.  le,  which  is  very 
similar  in  shape  to  the  core-level  photoemission  line  for  graphite 
(38).  Apparently  only  the  valence/core  excitonic  state  Is 
sufficiently  long  lived  to  significantly  participate  in  an  Auger 
process.  It  is  important  to  realize,  however,  that  the  lifetime 
effect  is  not  important  in  the  final  state,  and  all  of  the  dynamic 
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screening  excitations  could  possibly  produce  Auger  Intensity  In  the 
low-energy  region  of  the  llneshape. 

* 

The  position  of  the  structure  centered  at  about  240  eV  In  Fig. 
8  is  very  r’',r.r  that  expected  to  result  from  the  27  eV  plasmon 
sal  el’  tte  of  the  Auger  line  Itself.  Indeed,  much  of  the  Intensity 
In  this  region  (the  extrinsic  portion)  In  the  raw  data  of  Fig.  la 
wa3  removed  by  the  lo3s-deconvolution  procedure.  However,  in  the 
intrinsic  plasmon  emission  process  the  satellite  line  would  be 
expected  to  be  broader  than  the  Auger  line  itself,  i.e.,  it  would 
be  characterized  by  the  convolution  of  loss  structure  of  Fig.  1e 
with  the  true  Auger  llneshape  (72,73). 

The  width  of  the  240  eV  feature  can  be  estimated  by  assuming 
that  the  true  Auger  llneshape  is  well  represented  by  the  model 
spectrum  of  Fig.  8  and  subtracting  it  from  the  experimental 
llneshape.  The  result  of  this  procedure  13  shown  in  Fig.  9  where 
it  1 3  apparent  that  the  width  of  the  240  eV  feature  i3  much 
narrower  than  that  of  the  Auger  line.  Thus,  it  is  unlikely  that 
this  reature  can  be  identified  as  a  normal  intrinsic  plasmon 
satellite,  i.e.,  as  a  result  of  core-hole  filing. 

An  alternate  interpretation  is  suggested  by  the  recent  work  of 
Cini  and  D' Andrea  on  the  dynamical  aspects  of  screening  in 
determining  Auger  llneshapes  in  solids  (64).  In  this  work, 
dynamical  effects  (restricted  to  plasmon  emission)  are  included  by 
considering  a  plasmon  field  coupled  to  the  two  final-state  holes 
Interacting  through  the  bare  (unscreened)  hole-hole  repulsion, 
i.e.,  the  Inclusion  of  plasmon  emission  intrinsic  to  the  two-hole 
final  state  rather  than  the  core  hole.  The  physical  picture  of  the 
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process  Is  that  the  Auger  final-state  holes  can  be  created  with 
essentially  a  bare  hole-hole  repulsion  and  subsequently  become 
delocalized  through  the  emission  of  a  plasmon. 

Cini  and  D’Andrea  {61)}  find  that  this  new  channel  for  the 
decay  of  the  localized  holes  gives  rise  to  several  effects:  (1)  it 
can  result  in  plasmon  satellites  in  a  case  like  graphite  where  the 
spectrum  is  distorted  but  still  "bandlike",  (2)  it  can  produce 
broadened  but  characteristic  features  within  the  "one-electron" 
lineshape  if  the  hole-hole  distortion  is  strong  but  not  discrete 
and  (3)  it  can  produce  a  broadening  of  discrete  or  atomic-like 
features.  In  addition,  it  was  found  for  case  (l)  that  the  shape  of 
the  plasmon  satellite  did  not  reflect  the  shape  of  the  principal 
Auger  structure.  In  fact,  for  an  assumed  square  model  DOS,  Cini 
and  D'Andrea  !6if)  find  a  plasmon  satellite  that  is  a  broad  shoulder 
reminiscent  of  the  2^0  eV  feature  in  our  graphite  results  shown  in 
Fig.  2.  This  type  of  phenomenon  seems  a  likely  candidate  to 
account  for  the  remaining  discrepancies  between  our  model  and  the 
experimental  data  shown  in  Fig.  8. 

V.  SUMMARY  AND  CONCLUSIONS 

We  have  obtained  the  Auger  spectrum  for  P0C0  graphite  from 
data  taken  in  two  separate  laboratories  and  on  three  different 
types  of  el ect ron- energy  analyzers.  These  separate  data  were 
independently  corrected  for  the  effects  of  secondary-electron 
background  and  extrinsic  losses.  In  addition,  considerable  care 
was  taken  to  ensure  that  the  absolute  energy  scale  was  accurate. 

The  result  is  a  graphite  Auger  lineshape  which  we  feel  is  as  free 
as  possible  from  experimental  and  data-reduction  artifacts. 
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In  characterizing  the  Auger  structure  In  terms  of  the  one- 
electron  approximation  with  atomic  Auger  matrix  elements,  we  find  a 
model  lineshape  which  differs  considerably  from  that  of  the 
experiment.  Intensity  is  missing  In  the  model  function  at  both  the 
high-  and  low-energy  ends  of  the  spectrum  and  there  are  significant 
differences  in  relative  intensity  throughout  the  main  body  of  the 
line. 

To  characterize  the  discrepancies  seen  between  the 
experimental  lineshape  and  the  one-electron  Auger  model,  we  have 
considered  both  the  static  and  dynamic  aspects  of  inltlal-and 
final-state  screening.  We  find  that  the  static  polarization  effect 
of  initial-state  screening  has  a  negligible  influence  on  the 
lineshape.  However,  valence  electron  shake  up  into  the  core- 
excitonic  level  places  charge  in  an  energy  region  where  very  little 
exists  in  ground-state  graphite  giving  rise  to  significant  new 
Intensity  in  the  Auger  lineshape  just  below  the  Fermi  level. 
Modeling  the  dynamic  Auger  effect  by  the  inclusion  of  a  delta- 
function  density  of  states  at  the  Fermi  energy  and  assuming  that 
the  valence/ core  excitonic  electron  participates  in  the  Auger 
process  along  with  a  valence  electron,  results  in  a  dramatic 
improvement  between  the  measured  and  model  lineshapes  in  the  high- 
energy  region  for  an  effective  occupancy  in  this  excited  state  of 
0.27. 

The  distorting  effect  on  the  predicted  lineshape  resulting 
from  the  hole-hole  interaction  in  the  Auger  final  state  has  been 
modeled  using  the  Cini  expression  (35,36).  We  have  assumed  that 
the  empty  portions  of  the  o  and  v  band3  are  separated  sufficiently 
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from  each  other  and  from  the  filled  portions  to  permit  the  use  of 
the  Cini  filled-band  formalism  by  including  screened  hole-hole 
repulsion  parameters  for  the  o*o,  o*u^and  n*ir  contributions.  Under 
these  assumptions,  the  application  of  the  Cini  expression  results 
in  considerable  improvement  of  the  model  llneshape  in  the  region 
below  the  principal  maximum. 

The  final  area  of  disagreement  in  the  model  llneshape  consists 
of  a  shoul der- 1 1 ke  feature  on  the  low-energy  side  of  the  Auger  line 
which  is  not  accounted  for  by  localization  effects.  We  suggest,  on 
the  basis  of  more  recent  work  by  Cini  { 6 ti } ,  that  this  structure  is 
due  to  a  plasmon  effect  intrinsic  to  the  two-hole  final  state  in 
the  Auger  process.  The  adequate  characterization  of  this  feature 
will  undoubtedly  require  a  theoretical  model  which  includes 
multiple,  partially-filled  bands  with  the  inclusion  of  dynamical 
final-state  effects. 
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FIGURE  CAPTIONS 

(a)  The  raw  data  for  the  the  electron-excited  C(KVV) 

Auger  region  of  POCO  graphite  taken  in  the  N(E)  mode,  (b) 
The  raw  Auger  data  corrected  for  a  linear  secondary- 
electron  background,  (e)  The  electron  backscatter  spectrum 
of  POCO  graphite  using  290  eV  incident  electrons.  Curves 
(d)  and  (c)  show  a  comparison  of  the  low-energy  loss 
structure  of  (e),  on  an  expanded  vertical  scale,  and  the 
background  corrected  Auger  data  of  (b),  respectively. 

The  C(KVV)  lineshape  for  POCO  graphite  corrected  for 
secondary-electron  background  and  the  effects  of  extrinsic 
loss  processes  by  loss  deconvolution.  The  Fermi  level 
(FL)  in  this  kinetic  energy  scale  i3  located  at  the  C ( 1 s ) 
binding  energy  of  284.6  eV. 

A  broad  scan  high  resolution  XPS  spectrum  for  POCO 
graphite  showing  the  valence  and  Auger  regions  as  well  as 
the  C(1s)  line.  The  insets  show  expanded  views  of  the 
valence  and  Auger  regions. 

(a)  The  empirical  graphite  partial  densi ty-of-states 
( P DOS )  components  obtained  following  the  procedure  of 
Murday,  et  al .  {21}.  (b)  The  two-band  PDOS  where  the  o 
band  is  formed  by  summing  the  o  and  o  PDOS.  (c)  The 

3  P 

total  DOS  formed  by  summing  the  PDOS  of  (a).  The  vertical 
line  marked  FL  shows  the  position  of  the  Fermi  level  and 
the  curves  have  all  been  shifted  upward  in  energy  by  the 
C(1s)  binding  energy  of  284.6  eV.  The  dotted  portion  of 
(c)  shows  the  truncated  total  DOS  as  given  by  Murday,  et 
al .  {21  } . 

A  comparison  of  the  experimental  Auger  lineshape  for  POCO 
graphite  (solid  curve)  with  the  one-electron  model 
calculated  as  the  self-convolution  of  the  PDOS  from  Fig. 
4a,  modulated  by  the  symmetry-dependent  Auger  matrix 
elements.  The  vertical  line  marked  FL  shows  the  position 
of  the  Fermi  level.  The  dashed  portion  of  the  model  curve 
shows  the  effect  of  using  the  truncated  total  DOS  of 
Murday,  et  al .  {21}  from  Fig.  4c. 

A  comparison  of  the  effect  on  the  Auger  model  of  static 
initial-state  screening  using  the  partial  band  occupancy 
values  of  Dunlap,  et  al .  (solid  curve)  {26},  which  assumes 
that  only  the  it  band  is  involved  in  the  screening,  and  the 
occupancy  values  of  Binkley  {49}  which  includes  the 
contribution  from  all  bands  (dotted  curve).  The 
experimental  spectrum  from  Fig.  2  is  shown  as  the  dashed 
curve . 


Fig.  7.  A  comparison  of  the  experimental  Auger  lineshape  of  POCO 
graphite  (solid)  with  a  model  which  includes  the  initial- 
state  occupancy  values  of  Binkley  ( H 9 }  and  the  effect  of 
the  valence/core  excitonic  state  located  at  the  Fermi 
level  effectively  containing  0.27  electrons  (dashed). 

Fig.  8.  A  comparison  of  the  experimental  Auger  lineshape  of  POCO 
graphite  (solid  curve)  with  a  model  (dashed  curve)  which 
includes  initial-state  screening  (with  occupancy  values 
calculated  by  Binkley  (49)),  the  Fermi-level  valence/ core 
excitonic  state  and  the  hole-hole  interaction  distortion 
‘'rough  the  use  of  the  Cini  expression  (35). 
The  components  A,  B,  C  and  D  which  sum  to  the  model 
spectrum  are,  respectively,  the  o*o,  o*n,  n*n  .,j 
( o+x) *valence/core  exciton  {56}  contributions. 

Fig.  9.  The  difference  spectrum  resulting  from  the  subtraction  of 
the  model  spectrum  from  the  experimental  Auger  lineshape 
of  Fig.  8.  The  features  located  near  240  and  255  eV  both 
appear  considerably  narrower  than  the  Auger  lineshape 
itself. 
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